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1 . Introduction 

The  systems  of  distributions  (for  a random  variable  Y)  defined  by 


S : Z = y + 6 8jv  Y 

L 

(0  < Y) 

(1.1) 

SR:  Z = y ♦ 6 £n{Y/ (1-Y) } 

(0  < Y < 1) 

(1.2) 

S : Z = y + 5 sinh1  Y 

(1.3) 

with  6 > 0,  and  2 a unit  normal  variable,  have  been  described  by  Johnson 
(1949);  and  with  Z a standard  Laplace  (double  exponential)  variable,  (S£,  S^, 

Sy)  by  Johnson  (1954).  In  the  present  paper  we  study  analogous  systems, 
generated  by  ascribing  to  Z a standard  logistic  distribution  with  probability 
density  function  (pdf) 

fz(z)  = ez(l  ♦ ez)"2  (1.4) 


or,  equivalently,  with  cumulative  distribution  function  (cdf) 

Fz(z)  = (1  + e'2)"1  . 


(1.5) 


The  percentile  function  (inverse  of  cdf)  is 
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z = Hn{Fz(z)/(l  - F2(z))}.  (1.6) 

We  will  denote  these  new  systems  of  distributions  by  L^,  L^,  according 
us  transformations  (1.1),  (1.2),  (1.3),  respectively,  are  used.  In  view  of  the 
ckseness  of  the  shapes  of  the  logistic  and  normal  distributions,  it  is  to  be 
expected  that  the  new  systems  will  exhibit  some  similarity  in  shape  to  S^,  Sg, 
and  Sy. 

in  fact  since  (Johnson  and  Kotz  (1970),  p.6) 

|{1  + expf-nx/v^)  f1  - <!>(—)  | < 0.01  for  all  x (1.7a) 

t 2 

where  <Hu)  = (/2tT)  1 fUm  e ^ dt  the  difference  between  the  cdf  of  g (J  with 
parameters  y,  6 and  1.^  g ^ (with  the  same  values  of  £ and  X)  with  parameters 
Cy , ip6 > where 


15 

16 


1.7 


(1.7b) 


cannot  exceed  0.01.  Although  this  gives  good  agreement  in  the  central  part  of 
the  distribution,  there  can  be  gross  disparities  in  percentile  points  in  the 
tails. 


If  an  S curve  and  an  L curve  are  fitted  using  the  same  first  four  moments, 
we  do  not  find  tlat  (i)  the  ratios  (L/S)  for  y values  and  for  6 values  will  each 
be  about  1.7,  and  (ii)  the  values  of  £ and  X will  be  the  same.  Section  5.4 
(Table  4)  contains  an  example  wherein  this  is  clearly  illustrated.  There  can  be 
very  considerable  disparity  between  the  percentiles  in  the  tails  of  the  distribu- 
tions even  though  (1.7a)  is  valid. 

The  very  simple  formula  (1.6)  for  z in  terms  of  Fz(z)  is  an  obvious  practical 
advantage  of  the  new  systems.  Percentile  points  of  fitted  distributions  can  be 
obtained  very  simply. 
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2.  The  Log-Logistic  (LL)  Distribution 

This  distribution  has  been  studied  by  Shah  and  Dave  (1963).  The  pdf  is 


fY(y)  = 6 e'Y  y6*1  (1  ♦ e"Y  y6)"2  (y  >.  0;  6 > 0) 


(2.1) 


from  which  it  can  be  seen  that  it  is  a special  case  of  Burr's  (1942)  Type  XII 
system  of  distributions.  Dubey  (1966)  called  it  the  Weibull-exponential  distri- 
bution, and  fitted  it  to  some  business  failure  data  presented  by  Lomax  (1954). 

The  pdf  is  unimodal.  The  mode  is  at  y * 0 for  6 < 1 (giving  * reversed 
J- shaped  curve);  it  is  at  y = e Y (6-1)/ (6+1)  if  6 > 1.  The  cdf  is 


and  the  inverse  cdf  is 


rY(y)  * 1 - (1  ♦ e“Y  yV1 


y = en  [Fy(y)/(l-FY(y))]1/6 


(2.2) 


(2.3) 


where  ft  = y/6  . 

From  (1.1)  and  (1.4),  the  r-th  moment  of  Y about  zero  is 
u;  = p_m  er(z-Y)/6  e2  (1  + eV2  dz  . 

Putting  u * ez  (1  + e2)-1  so  that  du/dz  = ez(l  + e2)-^  and  e2  = u/(l-u) 

, -rft  fl  r/6  ,,  ..-r/6  , -rft  --1  . .-K 

y^,  = e / o u (1-u)  du  = e B(l+r6  , l-r6  ) 

-pQ 

= e r9  cosec  r0 

with  0 = tt/6,  provided  r < 6.  (If  r >.  6,  y^,  is  infinite.) 

In  particular 


E [Y]  = e"“  0 cosec  0 


(2.4) 


(2.5) 
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and 

-20  2 

Var(Y)  = e 0 cosec  0 (tan0  - 0)  . (2.6) 

The  r-th  central  moment,  pr>  is  a multiple  (depending  on  0,  but  not  on  fl) 

of  exp(-rfi),  and  so  the  moment-ratios  p /p^2  - in  particular  = p^/p^^  a°d 

2 

~ ^4/^*2  " <*°  not  depend  on  ft,  but  only  on  0.  The  (/B^^)  points  lie  on  a 

line  which  passes  through  the  'logistic  point'  (0,4.2)  (corresponding  to 

0 -►  0 (5  > 00)).  Figure  1 shows  that  ('log-logistic')  line,  and  also  the 

(lognormal)  line. 

The  (/d^*,  32)  points  for  Lg(Ly)  distributions  lie  'above'  ('below')  the 
line  in  Figure  1.  This  is  analogous  to  the  situation  for  S , S and  the  lognormal 
line. 

Table  1 gives  values  of  for  selected  values  of  /Bj,  for  points  on  the 

line. 

3.  The  I.  System 

If  transformation  (1.2)  is  valid,  with  Z standard  logistic,  then  the  pdf  of 

Y is 

fy(y)  = (5  eY  y6'1  (1-y)5'1  {(1-y)6  + eY  y6}'2  (0  < y < 1)  . (3.1) 

If  6 = 1,  we  get  a 'power  function'  pdf 

fy(y)  = - y + eY  y)“2  (0  < y < 1)  . (3.2) 

If,  in  addition,  y = 0 we  have  a standard  uniform  distribution. 

In  contrast  to  the  situations  when  Z has  a normal  or  Laplace  distribution, 

Lg  curves  cannot  be  multimodal.  There  is  a single  mode  (if  5 > 1)  or  antimode 
(if  <5  < 1)  at  the  unique  value  of  y between  0 and  1 satisfying  the  equation 


5 


The  right  hand  side  of  (3.3)  is  equal  to  1 when  y - whatever  the  value  of  6. 

The  derivative  of  the  logarithm  of  the  right  hand  side  is 


6(1- 

62)  y' 

-1,,  .-1 
U-y) 

{62  - (l-2y)2} 

• 

(3 

If  5 > 

1,  this 

is  negative, and  so 

the  mode 

is  at  y < % according  as  y < 

0. 

' 

For  6 < 1 

, the  curve  is 

U-shaped,  and 

the  antimode  is  between  ^(1-6] 

and 

*s(l+S). 

Since 

(3.4)  is 

positive  for  ^(1-6) 

< y < *s(  1-*- 6) 

it  follows  that 

the 

! 

i 

antimode  is  at  ^(1-5)  < 

y < 

\ for 

y < 0 and 

H < y < Js(i+6) 

for  y 

> 0. 

i 

1 

t 

TABLE  1: 

The  Log-logistic  Line 

/&T 

e2 

6 

'HI 

62 

6 

02 

6 

0.05 

4.21 

174.2 

0.55 

4.96 

16.26 

i.i 

7.47 

8.73 

0.10 

4.22 

87.14 

0.60 

5.11 

14.98 

1.2 

8.16 

8.13 

0.15 

4.26 

58.16 

0.65 

5.27 

13.90 

1.3 

8.94 

7.63 

0.20 

4.30 

43.69 

0.70 

5.45 

12.98 

1.4 

9.81 

7.21 

0.25 

4.35 

35.02 

0.75 

5.65 

12.18 

1.5 

10.79 

6 . 86 

0.30 

4.42 

29.26 

0.80 

5.86 

11.49 

1.6 

11.88 

6.55 

0.35 

4.50 

25.15 

0.85 

6.08 

10.88 

1.7 

13.11 

6.28 

0.40 

4.60 

22.08 

0.90 

6.32 

10.35 

1.8 

14.47 

6.04 

0.45 

4.71 

19.70 

0.95 

6.58 

9.87 

1.9 

16.00 

5.84 

0.50 

4.83 

17.80 

1.00 

6.86 

9.45 

2.0 

17.71 

5.65 

- 

(For  ✓'6^'  = 0, 

02  = 4.2 

and  6 

= "0 

. 


H 
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The  f/FTj",  R.,)  line  corresponding  to  6 = 1,  which  is  the  'lower'  boundary  of  the 
region  with  U-shaped  curves,  is  shown  in  Figure  1.  Apart  from  (3.2),  there  are 
no  J - or  reverse- J shaped  curves  in  the  L_  system. 

Following  an  analysis  similar  to  that  leading  to  (2.4)  we  obtain 

u'  = T U ♦ e-(z-^)/6rr  eZ(l  ♦ ezr2  dz 

r * — oo 

= /J  (1  ♦ e12  (1  - u)1/<S  u'1/6}'r  du  (3.S) 

with  ft  = y/6,  as  in  Section  2. 

Generally,  the  integral  in  (3.S)  must  be  evaluated  by  quadrature.  (Of 

course,  for  some  special  cases  - such  as  6 = 1 - explicit  solutions  can  be 

*11/6 

obtained.)  Expansion  of  the  integrand  as  a power  series  in  (1-u  ) ' cannot 

be  valid  over  the  whole  range  of  integration.  Dichotomy  of  the  interval  (0,1) 
according  as  (1-u-2)1^  < e$i  (i.e.  u % (1+e*^)’1)  leads  to  valid  expansions 
but  the  resulting  incomplete  beta  functions  must  themselves  (in  general)  be 
evaluated  by  quadrature. 


4.  The  Ly  System 

If  the  transformation  (1.3)  is  valid,  the  pdf  of  Y is 

f fy)  = AJLY_.  . W(y2+1))6 

Y (y2+l)'i  (l+eY(y+  /(y2+l))^}2 


(4.1) 


This  curve  is  unimodal,  with  mode  at  the  unique  value  of  y satisfying  the  equation 


Ml  - eY(y  + /(yZ+l)}  = y(y2+l)"'2 


(4.2) 


Note  that,  since  6 > 0,  the  left  hand  side  of  (4.2)  is  a monotonic  decreasing 
function  of  y;  the  right  hand  side  is  monotonic  increasing. 
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The  r-th  moment  of  Y about  zero  is 


K = 2'r  /T.  {eU'Y)/6  - e-(z^)/6}r  eZ(l  ♦ eZ)’2  dz 

= 2-r  l (-l)j  (r)  e" (r_2jKl  £ e(r-2j)z/6  ^ + ^-2  dz 

j=0  J 

r 

= 2"r  | (-1)-*  (?)  e"^r-2^fi  (r-2j) 0 cosec(r-2j)0 

j-0  J 


(4.3) 


provided  r < 6.  If  r > 6,  is  infinite.  When  r = 2j , '0  cosec  O'  is  inter- 

preted as  lim  6 cosec6  = 1.  (As  in  Section  2,  ft  = y/6  and  6 = tt/6.) 

0-*O 

Since  (-a)cosec(-a)  = a cosec  a,  (4.3)  can  be  simplified. 


For  r even: 


!ir- 1 

u'  = 2 r(-l)  iT  (Jj  + 2 ^ 1 (-1)^  (T)(r-2j)0  cosec(r-2j)0  cosh(r-2j)ft 

^ j=0  J 


For  r odd: 


In  paiticular 


-fr-n  i r 

U'  = -2  1 l (-1)J  (.)(r-2j)0  cosec(r-2j)0  sinh(r-2j)ft 

j=0  3 


pj  = -0  cosec  0 sinh  ft 


p£  = 0 cosec  20  cosh  2ft  - j 


p^  = - ^ 0(cosec  30  sinh  3ft  -cosec  0 sinh  ft) 

pj  = -J  0(cosec  40  cosh  4ft  - 2 cosec  20  cosh  2ft  ) ♦ 


(4.4) 


(4.5) 


(4.6) 


so  that 


E[Y]  = -0  cosec  0 sinh  ft 

Var(Y)  = Jj{  (0  cosec  0)2-l}  + !50  cosec2  0(tan0  - 0)cosh  2ft 


(4.7) 


I ■!■■■ 


.....  - ■ i 
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5 . Fitting  the  Distributions 

5.1.  General . The  methods  of  fitting  described  by  Johnson  (1949)  for  the  S 
systems  are  also  applicable  to  the  L systems,  with  the  advantage  that  the  cdf 
of  the  logistic  is  simpler  than  that  of  the  normal  distribution. 

Introducing  the  location  and  scale  parameters  £,  X respectively,  by  the 
transformation  Y = (X-£)/A  we  obtain  a three  parameter  family  for  L^: 

Z = y + <5  £n  (X-O  (X  > £)  (5.1) 

and  four  parameter  families 

Z = Y + <5  in  (iTPx)  U < X < £ ♦ X)  (5.2) 

Z = y + 6 Jlinh"1  (^£)  (5.3) 

for  Lg,  Ly  respectively. 

tfe  consider  fitting  by  the  methods  of  moments,  percentile  points  and 
maximum  likelihood. 

Moments.  If  four  parameters  are  to  be  fitted,  the  first  four  moments  of 
the  distribution  of  X are  equated  to  those  of  the  fitted  curve.  It  is  convenient 
to  use  them  in  the  form  of  mean  (tip,  the  variance  (p^)  and  the  moment  x-atios 
(=  Uj/t-/  ) and  p.,  (=  p^/p") . Since  and  are  determined  by  y>  <5  and 
(for  Lg  and  L^,  as  for  Sg  and  S^)  conversely,  the  first  step  is  to  determine 
Y,  from  the  specified  values  of  and  82*  This  could  be  done  using  special 
tables  (Table  1,4),  or  of  course  a computer  program.  Once  these  values  are  found 
the  values  of  £ and  X are  determined  from  the  equations 


F-:[X|y,  6]  = p ♦ AE[Y|y.  6] 
var(X | y,  6)  = X2  var(Y|Y,  6) 


(5.4) 
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Percentile  Points.  If  four  parameters  are  to  be  fitted,  then  (using  (1.6)), 
four  equations  of  form 


P.  XP.~  ^ 


2 = ^n(r^-)  = y + 6 f(— i ) (j  = 1,2, 5,4) 


(5.5) 


where  Xp  is  the  estimated  100  P.%  point  of  the  distribution  of  X,  have  to  be 

j _ _ _ _ J 

solved  for  y,  <$,  £ and  A. 

Solutions  of  especially  simple  form  are  obtained  by  taking  symmetric  per- 
centiles . ith  P.  = 1-P_ ; P,  = 1-P,  . 

3 2'  4 1 

It  is  interesting  to  note  that  with  this  choice  of  P's,  provided  the  distri- 
bution of  Z is  symmetrical  about  zero  then  the  solutions  of  (5.5)  for  f,  and  A 
must  be  linked  by  the  same  equation 


XP  Xl-P"?  XP 

f(-4— ) + f( — ) = f(— r 


Xl-P,-^ 

) * f( r~ — ) 


(5.6) 


whatever  the  actual  distribution  of  Z. 
In  fact,  from 


6 f(— 4— ) 

A 


Xl-P.^ 


-Zl-P.  = "Y 
J 


- 6 f( J— ) (j  = 1,2) 

A 


it  follows  that  the  common  value  of  each  side  of  (5.6)  is  -2Q  . 
Another  relation  is 

Xl-P  XP?"^  Xl-P  XP 

z (f(-----)  - f(— =— = z.  p (fC — ) - f(-4— )1  • 

1 A A 2 A A 


(5.7) 


This  relation  does  depend  on  the  actual  distribution  of  Z.  Equations  (5.6) 
and  (5.7)  suffice  to  determine  £ and  A. 

In  certain  cases  (e.g.  for  L„  - see  below)  equation  (5.6)  can  be  solved 
explicitly  giving  A as  a function  of 
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If  £ and  X are  known,  two  percentile  points  suffice  to  determine  y and  6. 
Knowledge  of  t,  and/or  X is  more  usual  wnen  a hounded  curve  is  being  fitted. 

Fitting  of  ^ in  these  cases  is  discussed  in  Section  5.3. 

Max; '..am  Likelihood.  If  E,  and  X are  known,  the  transformed  values  of 
f((x-£,)/X)  can  be  used,  and  methods  appropriate  to  fitting  logistic  distributions 
(e.g.  . Johnson  and  Kotz  (1970,  Chapter  23))  can  be  applied.  If  £ and  X are  not 
known,  a series  of  pairs  of  values  of  these  parameters  can  he  used  and  the 
corresponding  maximized  likelihoods  calculated,  the  maximum  likelihood  estimates 
being  found  by  trial  and  error  - or,  of  course,  a computer  program  which  directly 
maximizes  the  likelihood  function  can  be  used.  Except  in  special  cases,  maximum 
likelihood  calculations  are  rather  lengthy,  but  not  impracticably  so. 

5.2.  Fitting  the  System.  If  £ is  known,  the  transformed  values 

Y.  = log(X.-f,) 

can  be  used,  and  the  parameters  y and  6 fitted  by  methods  appropriate  to  fitting 
a logistic  distribution. 

When  £ is  not  known  there  are  three  parameter  values  to  be  fitted  (£*Y>  and  6) . 
With  the  method  of  percentile  points,  it  is  natural  to  use  the  sample  median 
Xp  j.  and  the  upper  and  lower  sample  100  P%  points,  p and  Xp  respectively. 

Noting  (1.6),  this  leads  to  the  equations 

0 = y + 6 MX0  5 - l) 

-«.»(!>/ (1-1’)}  - Y ♦ <5  «,n(Xl  p - Q (5.8) 

Zn{P/(l-P)}=  Y + 6 ?,n  ( Xp  - l) 
which  have  the  unique  solutions 


i 
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£ ■ 


X0.5  " Vp^P 


2X0.S  “ Xl-P  ' XP 


-2  £,n{P/  (1-P) } 
an{(Xl  p - 0(Xp  - £)> 

Y = -6  £n(X0>5  - l) 


A 

provided  £ < Xj,. 

Of  course, several  values  of  P may  be  used.  The  resulting  estimates  Y 
and  6 should  be  reasonably  consistent  with  each  other  if  an  distribution  is 
suitable.  If  only  one  set  is  to  be  used,  P = 0.10  (i.e.  the  sample  median, 
anJ  upper  and  lower  sample  deciles)  would  seem  to  be  a good  choice. 

3/2 

We  now  come  to  the  method  of  moments . The  skewness  parameter  /£j  = 


is 


/b^ToT 


2 3 3 

36  cosec  36  - 66  cosec  6 cosec  26  + 39  cosec  6 

{9  cosec20  (tan6  - 6)  }^2 

2 2 3 

_3  cosec  39-9  cosec  9 sec  9 + 9 cosec  9 
cosec^9  (tan9  - 9)^2 


(5.10) 


3/2 

Equating  (5.3)  to  v^"  = m^/i^'  » the  samPle  skewness  ratio,  and  solving  for 
6 gives  the  moment  estimator,  9,  of  9 . The  function  (5.3)  is  an  increasing  func- 
tion of  9,  and  tends  to  infinity  as  9 + ® , so  provided  > 0,  there  will  be  a 
unique  value  for  6.  The  estimator  for  <5  is 


6 = tt/9  . 


Example  1.  As  can  be  seen  from  Figure  1,  the  line  is  in  the  Pearson  Type  IV 
region.  Suppose  we  want  to  use  an  curve  to  approximate  a Type  IV  with 
= 0.4,  62  = 4.6.  From  Table  l (with  extra  decimal  places),  6 = 22.0803 
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: 


yields  /Bj  = 0.4,  B2  4.5991  with  E[Y|y=0,  6)  = 1.0034,  6(Y |y=0,  6)  = 0.0828 
whence  y = -55.0,  f,  = 12.1184.  Standardized  percentiles  of  the  fitted  distri- 
bute are  compared  with  those  of  the  Type  IV  distribution  (obtained  by  inter- 
polation in  the  tables  of  Bouver  and  Bargmann  (1974))  in  Table  2. 

Table  2.  Comparison  of  Standardized  Type  IV  and  Percentiles  (Bj  = 0.4,  = 4-6) 


% 

0.1 

0.25 

0.5 

1 

2.5 

5 

10 

25 

50 

Type 

IV 

-3.315 

-2.887 

-2.574 

.2.262 

-1.854 

-1.532 

-1.186 

-0.641 

-0.046 

ll 

-3.285 

-2.910 

-2.615 

-2.310 

-1.887 

-1.548 

-1.185 

-.627 

-.041 

% 

75 

90 

95 

97.5 

99 

99.5 

99.75 

99.9 

Type 

IV 

0.587 

1.237 

1.686 

2.127 

2.717 

3.182 

3.674 

4.346 

0.575 

1.223 

1.682 

2.139 

2.753 

3.231 

3.723 

4.395 

5.3.  Fitting  the  System. 

Moments.  If  all  four  parameters  (£,A,y,6)  have  to  be  estimated,  it  is  neces- 
sary to  have  special  tables  (i)  to  determine  y,  6 from  /jTj’  and  B2  an<*  (ii)  giving 
values  of  E[Y|y,  6),  and  S.D(Y|y,  6)  for  these  values  of  y and  6.  Such  tables 
are  in  preparation  but  as  yet  they  are  not  sufficiently  extensive  for  practical  use. 

Percentile  Points.  For  Lg  curves,  introducing  the  notation  £p  = Xp  - 


equations  (5.6)  and  (5.7)  become 
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respectively.  From  (5.11) 


,\  = 


Cp  p (Sp  ♦ S,  p ) • Kp  Kt  p (Cp  ♦ Si  p ) 
F1  1_Fi  l2  1_P2  P2  l'y2  *1  l'yl 


(5.13) 


^1-P 
*1  1 yl 


tp  h-p 
l2  1 *2 


These  equations  may  be  solved  by  a process  of  trial  and  error,  illustrated 
in  Example  2(a)  below. 

Techniques  of  solving  similar  equations  for  fitting  curves,  devised  by 

Bukac  (1972),  Mage  (1978)  and  Shapiro  (1978),  may  also  be  used  for  Lg.  Since 

Pj  and  ?2  may  be  chosen  arbitrarily,  it  is  possible  to  arrange  that  Zj_p  /Zj_p 

is  an  integer.  We  do  this  in  Example  2(a),  though  it  is  not  necessary  to  do  so. 

v 

Devices  for  simplifying  estimation  of  the  four  parameters  for  S^  (Bukac 
(1972),  Mage  (1978),  Slifker  and  Shapiro  (1979))  can  equally  be  applied  to 
fitting  L_.  We  do  not  pursue  these  matters  here. 

D 

In  Example  2(b)  and  2(c)  we  fit  the  same  data  as  in  2(a),  when  both  £ and 
X are  assumed  known,  and  when  £ (but  not  X)  is  assumed  known. 

If  the  range  of  variation  (and  so  both  £ and  X)  be  known,  then  only  two  sample 
percentiles  are  needed  to  estimate  y and  6.  It  is  convenient  to  choose  them 
symmetrically  and  use  the  equations 

- - 


•zi-p  ■ * * 6 to  (ft> 


(5.14) 


1-P 


:,-p  = " + 6 ^ (1R^) 


where  £p  = Xp  - whence 

6 = 

and 


2 Z 


1-P 


p U'Sp)  -1 

q.o-qjv  1 


(5. IS) 


Y = - 


1 ? „ *P  ^1-P 

2 6 (X-Cp)(X-C1_p) 


(5.16) 


- 
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If  only  one  end  point  - say  the  lower  (£)  - is  known,  then  three  sample 
percentiles  are  needed,  from  which  to  estimate  X,  y and  6.  It  is  convenient  to 
take  the  sample  median  XQ  in  addition  to  Xp  and  X^  p.  From  the  resulting 
equations  we  obtain 

( ^0-5  }2  = _ CP  Cl-P 
*"^o.5  (X-£p) (X-£lp) 

whence 


~ ^0.5  ^0.5  ^P  + ^l-P^  2 ^1-P^ 

* 2 

^0.5  ' ^P  ^1-P 

y and  6 are  then  calculated  from  (5.15)  and  (5.16). 


(5.17) 


Example  2.  (a)  We  will  use  the  data  in  Table  6-4  of  Hahn  and  Shapiro  (1967). 


For  fitting  four  parameters  we  use,  in  addition  to  the  values  X = 0.84, 

X = 1.42  employed  by  Hahn  and  Shapiro  for  fitting  S_,  the  values  X = 0.97 


‘.91 

A 

X ,0,0  = 1.18.  (These  values  are  chosen  to  make 
. bo  jo 

„ 0.91  _ „ 0.6838  - , . 

= "n  0.09  ~ 2-3136  ~ 3 n 0.3162  " 3 2 .6838 


‘.3162 


^.91 


The  trial  and  error  calculations  are  set  out  briefly  below. 


(1) 

(2) 

£PU-£i  P) 

£ 

HO 

g(0.09) 

g(0 . 3162) 

(D/(2) 

g(P)  = £n 

i i t 

0.5 

6.32 

-1.097 

-0.424 

2.59 

£i.pCX-£p) 

0.75 

0.95 

-3.130 

-1.009 

3.10 

We  want  g(0.09)/g(0. 3162)  = 3 . 

0.72 

1.12 

-2.631 

-0.886 

2.97 

0.73 

1.06 

-2.779 

-0.924 

3.01 

Hence 

we  obtain 

l = 0.73, 

X - 1.06. 

These  lead 

to  estimates 

y = 1.276,  6 

= 1.665 

Percentage  distribution  is  shown  in  column  (3)  of  Table  3. 

(b)  If  we  take  £ = 0.5,  X = 1.5  and  use  X^  ^ and  9.^  we  obtain 

y = 1.049,  6 = 2.740  . 
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Since  we  are  using  the  same  values  for  £ and  X as  Hahn  and  Shapiro  used  in 
fitting  an  Sg  curve  to  the  same  data  (and  so  the  same  Xg  and  Xg  gj),  the 
values  are  each  in  the  same  ratio  (2.3136/1.3408  » 1.726  - note  that 
4>(1.3408j  = 0.91)  to  the  corresponding  values  for  their  fitted  Sg. 

Percentage  distributions  of  the  fitted  are  shown  in  column  (S)  of  Table 
3.  Column  (4)  contains  values  obtained  by  Hahn  and  Shapiro  for  their  fitted  Sg. 

(c)  If  Xp  g = 1.07  is  used  in  addition  to  and  Xg  we  obtain  from 

(5.17) 

A = 4.36 

which,  of  course,  agrees  with  the  value  obtained  by  Hahn  and  Shapiro  when  fitting 
an  Sg.  The  values  for  y and  6, 

y = 7.614  and  6 = 4.019 

are  again  in  the  same  ratio  (1.726)  to  Hahn  and  Shapiro's  Sg  values. 

Percentage  distributions  of  the  fitted  L0  are  shown  in  column  (7)  of  Table 
3.  Column  (6)  contains  values  obtained  by  Hahn  and  Shapiro  for  their  fitted  Sg. 

Comparing  the  three  fits  (columns  (3),  (5)  and  (7))  it  is  noteworthy  how 
similar  shaped  curves  can  be  obtained  with  quite  different  values  for  y and  6, 
by  appropriate  adjustment  of  £ and  X (and  conversely) . 


S 
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Table 

3 - 

Comparison  of  Fitted 

SB  ***  LB 

Curves 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Production 

Time 

(.min) 

Actual 

% 

lb 

SB 

lb 

SB 

lb 

< 0.695 

0.9 

- 

0.9 

1.5 

0.4 

0.9 

0.695-0.795 

3.7 

3.7 

4.7 

4.2 

4.5 

4.2 

0.795-0.895 

12.6 

14.0 

10.3 

8.8 

12.8 

11.0 

0.895-0.995 

18.4 

18.8 

15.2 

14.6 

18.7 

18.4 

0.995-1.095 

18.8 

18.6 

17.4 

18.4 

19.1 

20.6 

1 .095-1.195 

15.8 

15.3 

16.9 

18.1 

15.7 

16.8 

1. 195-1 .295 

12.2 

11.3 

13.9 

14.3 

11.3 

11.2 

1.295-1.395 

7.6 

7.8 

10.2 

9.4 

7.4 

6.8 

1.39S- 1.495 

5.0 

5.1 

6.1 

5.5 

4.5 

4.0 

1.495-1.595 

2.8 

3.1 

3.1 

2.9 

2.6 

2.4 

1.595-1.695 

1.1 

1.7 

1.0 

1.4 

1.4 

1.4 

1.695-1.795 

0.9 

0.6 

0.3 

0.6 

0.8 

0.9 

a 1.795 

0.2 

- 

0.0 

0.2 

0.8 

1.5 

£ 

0.73* 

0.5 

0. 

,5 

X 

1.06* 

1.5 

4. 

,36* 

Y 

1.276* 

0.608* 

1.049* 

4.411* 

7.614* 

6 

1.665* 

1.587* 

2.740* 

2.328* 

4.019* 

"estimated 
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5.4.  Fitting  Ly  Distributions 

Moments:  Table  4 gives  values  of  5 and  ft(*  y/6),  corresponding  to 
^7  3 0.00(0.05)1 (. 1)2.0  combined  (for  each  /gj)  with  twenty  values  of  f?2 
increasing  by  intervals  of  0.2,  starting  from  a value  just  'below'  the  L,  line. 
The  simultaneous  nonlinear  equations 


(y,  6)  = 

62(y,6)  = 02 


(5.18) 


were  solved  for  y and  6 using  Brown's  (1967)  algorithm.  The  search  for  values  of 
n and  6 continued  until  | ^ (y, 6)  - /fj^l  < 10~6  and  |&2(y,6)  - 62|  < 10"6.  The 
use  of  the  table  is  illustrated  by  the  following  example. 


Example  3.  We  will  compare  a standardized  Pearson  Type  IV  curve  with  * 0.9, 
82  = 8.6  against  an  Ly  curve  having  the  same  first  four  moments. 

From  Table  4 we  obtain  6 « 6.0151  and  ft  * -O.S2SO,  where 
y = 6.0151  x (-0.5250)  * -3.1580.  Using  (4.7)  we  obtain 

E[Y]  = 0.5753;  S.D. [Y]  * 0.3712  . 


Using  (5.4)  (with  'observed'  mean  zero,  and  standard  deviation  1)  we  obtain 
X * (0.3712)"1  = 2.6940,  £ = 0 - (2.6940  x 0.5753)  * -1.5498.  So  the  fitted 

Ly  curve  is  defined  by  the  statement  that 

Z = -3.1580  ♦ 6.0151  sinh’1  {(X  ♦ 1.5498/2.6940} 
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Table  4:  To  Facilitate  Fitting  Curves 


- 0.00 


«2 

< 

Q 

V 

0 

4.40 

16.1153 

0.0000 

0.0000 

0.1136 

4.60 

11.7442 

0.0000 

0.0000 

0.1571 

4.80 

9.8648 

0.0000 

0.0000 

0.1884 

5.00 

8.7752 

0.0000 

0.0000 

0.2131 

5.20 

8.0510 

0.0000 

0.0000 

0.2337 

5.40 

7.5295 

0.0000 

0.0000 

0.2512 

5.60 

7.1336 

0.0000 

0.0000 

0.266S 

5.80 

6.8216 

0.0000 

0.0000 

0.2800 

6.00 

6.5687 

0.0000 

0.0000 

0.2920 

6.20 

6.3590 

0.0000 

0.0000 

0.3028 

6.40 

6.1822 

0.0000 

0.0000 

0.3126 

6.60 

6.0308 

0.0000 

0.0000 

0.3215 

6.80 

5.8998 

0.0000 

0.0000 

0.3297 

7.00 

5.7850 

0.0000 

0.0000 

0.3373 

7.20 

5.6837 

0.0000 

0.0000 

0.3442 

7.40 

5.5935 

0.0000 

0.0000 

0.3507 

7.60 

5.5127 

0.0000 

0.0000 

0.3567 

7.80 

5.4398 

0.0000 

0.0000 

0.3623 

8.00 

5.3738 

0.0000 

0.0000 

0.3676 

8.20 

5.3137 

0.0000 

0.0000 

0.3725 

- o.io 

®2 

6 

n 

u 

a 

4.40 

16.8772 

-0.1905 

0.1927 

0.1103 

4.60 

12.0026 

-0.1306 

0.1325 

0.1549 

4.80 

10.0035 

-0.1057 

0.1076 

0.1867 

5.00 

8.8648 

-0.0912 

0.0933 

0.2118 

5.20 

8.1148 

-0.0815 

0.0837 

0.2325 

5.40 

7.5779 

-0.0745 

0.0767 

0.2502 

5.60 

7.1720 

-0.0690 

0.0713 

0.2656 

5.80 

6.8530 

-0.0647 

0.0671 

0.2792 

6.00 

6.5950 

-0.0612 

0.0636 

0.2913 

6.20 

6.3815 

-0.0582 

0.0607 

0.3021 

6.40 

6.2017 

-0.0557 

0.0582 

0.3120 

6.60 

6.0481 

-0.0535 

0.0560 

0.3210 

6.80 

5.9150 

-0.0516 

0.0542 

0.3292 

7.00 

5.7987 

-0.0499 

0.0525 

0.3368 

7.20 

5.6961 

-0.0485 

0.0510 

0.3438 

7.40 

5.6048 

-0.0471 

0.0497 

0.3503 

7.60 

5.5230 

-0.0459 

0.0485 

0.3S64 

7.80 

5.4494 

-0.0448 

0.0474 

0.3620 

8.00 

5.3826 

-0.0438 

0.0464 

0.3673 

8.20 

5.3219 

-0.0429 

0.0455 

0.3722 

- 0.05 


*2 

6 

0 

V 

a 

4.40 

16.2968 

-0.0909 

0.0916 

0.1127 

4.60 

11.8071 

-0.0638 

0.0646 

0.1565 

4.80 

9.8988 

-0.0S20 

0.0529 

0.1879 

5.00 

8.7974 

-0.0451 

0.0461 

0.2128 

5.20 

8.0668 

-0.0404 

0.0414 

0.2334 

5.40 

7.5415 

-0.0369 

0.0380 

0.2510 

5.60 

7.1432 

-0.0343 

0.0354 

0.2663 

5.80 

6.8294 

-0.0322 

0.0333 

0.2798 

6.00 

6.5752 

-0.0304 

0.0316 

0.2918 

6.20 

6.3647 

-0.0290 

0.0302 

0.3026 

6.40 

6.1870 

-0.0277 

0.0289 

0.3124 

6.60 

6.0351 

-0.0266 

0.0279 

0.3214 

6.80 

5.9036 

-0.0257 

0.0270 

0.3296 

7.00 

5.7884 

-0.0249 

0.0262 

0.3372 

7.20 

5.6868 

-0.0241 

0.0254 

0.3441 

7.40 

S.S963 

-0.0235 

0.0248 

0.3S06 

7.60 

5.5153 

-0.0229 

0.0242 

0.3566 

7.80 

5.4422 

-0.0223 

0.0236 

0.3623 

8.00 

5.3760 

-0.0218 

0.0232 

0.3675 

8.20 

5.3157 

-0.0214 

0.0227 

0.3725 

*T" 

0.15 

*2 

6 

a 

V 

o 

4.40 

18.0083 

-0.3121 

0.3188 

0.1065 

4.60 

12.3522 

-0.2040 

0.2076 

0.1523 

4.80 

10.1867 

-0.1628 

0.1661 

0.1846 

5.00 

8.9813 

-0.1396 

0.1430 

0.2100 

5.20 

8.1972 

-0.1243 

0.1277 

0.2311 

5.40 

7.6401 

-0.1132 

0.1167 

0.2489 

5.60 

7.2211 

-0.1048 

0.1084 

0.264S 

5.80 

6.8931 

-0.0981 

0.1017 

0.2782 

6.00 

6.6285 

-0.0926 

0.0963 

0.2904 

6.20 

6.4102 

-0.0881 

0.0918 

0.3013 

6.40 

6.2266 

-0.0842 

0.0880 

0.3113 

6.60 

6.0698 

-0.0809 

0.0847 

0.3203 

6.80 

5.9344 

-0.0779 

0.0818 

0.3286 

7.00 

5.8161 

-0.0754 

0.0793 

0.3362 

7.20 

5.7118 

-0.0731 

0.0770 

0.3433 

7.40 

5.6191 

-0.0711 

0.0750 

0.3498 

7.60 

5.5361 

-0.0692 

0.0732 

0.3559 

7.80 

5.4614 

-0.0676 

0.0715 

0.3616 

6.00 

5.3938 

-0.0661 

0.0700 

0.3669 

8.20 

5.3322 

-0.0647 

0.0686 

0.3719 
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/T[  - 0.20 


» 2 

6 

0 

u 

a 

4.40 

20.0017 

-0.4873 

0.5089 

0.1019 

4.60 

12.9064 

-0.2896 

0.2966 

0.1406 

4.00 

10.4634 

-0.2259 

0.2313 

0.1817 

5.00 

9.1539 

-0.1917 

0.1967 

0.2077 

5.20 

0.3170 

-0.1696 

0.174S 

0.2290 

5.40 

7.7305 

-0.1539 

0.1S89 

0.2472 

5.60 

7.2920 

-0.1421 

0.1471 

0.2629 

5.00 

6.9509 

-0.1327 

0.1378 

0.2768 

6.00 

6.6767 

-0.1252 

0.1302 

0.2891 

6.20 

6.4512 

-0.1109 

0.1240 

0.3002 

6.40 

6.2621 

-0.1135 

0.1107 

0.3102 

6.60 

6.1009 

-0.1089 

0.1141 

0.3194 

6.00 

5.9620 

-0.1049 

0.1101 

0.3277 

7.00 

5.8408 

-0.1014 

0.1067 

0.3354 

7.20 

5.7341 

-0.0983 

0.1036 

0.3425 

7.40 

5.6393 

-0.0955 

0.1098 

0.3491 

7.60 

5.5546 

-0.0930 

0.0983 
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5.1969 

1.8380 

1 

18.20 

5.5494 

-2.0440 

4.0066 

1.4467 

1 

18.40 

5.5107 

-1.8815 

3.3860 

1.2471 

18.60 

5.4739 

-1.7629 

2.9905 

1.1228 

18.80 

5.4388 

-1.6702 

2.7107 

1.0367 

19.00 

5.4054 

-1.5944 

2.4997 

0.9731 

19.20 

5.3736 

-1.5308 

2.3333 

0.9240 

19.40 

5.3432 

-1.4760 

2.1979 

0.8849 

19.60 

5.3142 

-1.4282 

2.08SO 

0.8530 

19.80 

5.2864 

-1.3858 

1.9891 

0.8265 

20.00 

5.2597 

-1.3480 

1.9063 

0.8040 

20.20 

S.2342 

-1.3138 

1.8340 

0.7848 

20.40 

5.2097 

-1.2828 

1.7701 

0.7681 

20.60 

5.1862 

-1.2544 

1.7131 

0.7536 

20.80 

5.1636 

-1.2283 

1.6620 

0.7408 

21.00 

5.1418 

-1.2042 

1.6157 

0.7295 

, 

21.20 

5.1209 

-1.1818 

1.5736 

0.7193 

21.40 

5.1007 

-1.1609 

1.S3S1 

0.7103 

21.60 

5.0812 

-1.1415 

1.4998 

0.7021 

r 
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ha  ; a standard  logistic  distribution.  Percentile  points  of  the  Type  IV  curve 
(iron)  F.ouver  and  Bargmann  (1974))  and  this  Ly  curve  are  compared  in  Table  5, 
which  also  contains  values  for 

(i)  the  Sy  curve  with  the  same  firs*  four  moments,  defined  by 
Z*  = -0.4048  + 1.455  sinh'1  {(X  + 0 . 3842)/l .0765} 


having  a standard  normal  distribution,  and 

(ii)  the  Sy  curve  defined  by  1**  = Z/1.7  (see  (1.7  a,b)) 

We  note  that  (a)  the  Ly  curve  agrees  more  closely  with  the  Type  IV  than  does 
the  first  (moment- fi tted)  Sy  in  the  tails  and  (b)  although  the  probability  integral 
of  'he  ‘econd  Sy  curve  must  agree  with  the  Ly  to  within  0.01,  there  is  considerable 
divergence  between  the  percentiles  near  0 or  100. 


Table  5 . Comparison  of  Type  IV,  Ly  and  Sy  Standardized  Percentiles 


% 

0.1 

0.2S 

0.5 

1 

2.5 

5 

10 

25 

50 

Type  IV 

-3.345 

-2.841 

-2.492 

-2.164 

-1.748 

-1.437 

-1.114 

-0.619 

-0.078 

Lu 

-3.359 

-2.865 

-2.526 

-2.220 

-1.777 

-1.454 

-1.118 

-0.609 

-0.070 

Sy(i) 

-3.707 

-3.086 

-2.656 

-2.257 

-1.778 

-1.417 

-1.073 

-0.585 

-0.081 

vu> 

-2.S07 

-2.281 

-2.100 

-1.908 

-1.628 

-1.387 

-1.109 

-0.632 

-0.070 

75 

90 

95 

07.5 

99 

99.5 

99.75 

99.9 

Type  IV 

0.524 

1.192 

1.690 

2.211 

2.959 

3.586 

4.279 

5.316 

Lu 

0.520 

1.178 

1.677 

2.203 

2.967 

3.609 

4.318 

5.375 

Sy(i) 

0.490 

1.162 

1.686 

2.244 

3.051 

3.721 

4.451 

5.513 

Sy(ii) 

0.546 

1.166 

1.574 

1.955 

2.432 

2.781 

3.123 

3.571 

(i)  moment 

fit; 

(ii)  -1. 

8576  «■  3. 

.5383  sinh'1  {(X 

+ 1.5498)/2.6940}  unit 

normal 

i 
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